Mobility of Bloch Walls via the Collective Coordinate Method 



M.A. Desposito* 1 , A. Villares Ferrer 2 , A.O. Caldeira 2 



(^ 



X 



and A.H. Castro Neto' 



: J . 



q_) ■ 1 Departamento de Fisica, Facultad de Ciencias Exactas y Naturales, 

iy-s . Universidad de Buenos Aires, RA-1^28, Buenos Aires, Argentina. 

<n : 

Institute de Fisica "Gleb Wataghin" , 
q . Departamento de Fisica do Estado Solido e Ciencia dos Materiais, 

Universidade Estadual de Campinas, 13083-970, Campinas, SP, Brasil. 

3 Dept. of Physics, University of California, Riverside CA, 92521, USA. 
+-* 

a 

^ ■ Abstract 

a 
o 

q ' We have studied the problem of the dissipative motion of Bloch walls con- 

sidering a totally anisotropic one dimensional spin chain in the presence of 
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a magnetic field. Using the so-called "collective coordinate method" we con- 
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q^ ' excitations of the system. It allows us to analyze the Brownian motion of the 

+- > 
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struct an effective Hamiltonian for the Bloch wall coupled to the magnetic 



wall in terms of the reflection coefficient of the effective potential felt by the 
excitations due to the existence of the wall. We find that for finite values of 



h3 

the external field the wall mobility is also finite. The spectrum of the potential 

o" 

at large fields is investigated and the dependence of the damping constant on 



temperature is evaluated. As a result we find the temperature and magnetic 



C3 ■ field dependence of the wall mobility. 
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I. INTRODUCTION 

It is a well-known fact that ultimately due to magnetic dipole interaction, different do- 
mains are formed in magnetic systems 0. In many situations, the physical region separating 
two different magnetic domains — the domain wall — must be treated as a physical entity be- 
cause it has a characteristic behaviour when acted by external agents. For instance, it is 
known that the response of a magnetic system to a frequency dependent external mag- 
netic field depends on whether domain walls are present 0. Domain walls can also move 
throughout the system and this motion happens to be dissipative M. 

A particularly interesting kind of domain wall is commonly found in low dimensional 
ferromagnetic systems. These are the so-called Bloch walls @]. It is known that these walls 
perform dissipative motion || due to the presence of the elementary excitations which can 
be scattered by the wall as it moves and the momentum transferred to them reduces the 
speed of the wall. 

The primary aim of this work is to study the influence of finite temperatures in the 
mobility of these Bloch walls. For this purpose we start by considering a microscopic model 
for a one dimensional ferromagnet containing hard and easy-axis anisotropies and subject to 
an external magnetic field. A semiclassical picture provides us with the localized solutions 
for the spin configurations which are the solitons corresponding to the walls. 

Making use of a recently developed method for the analysis of the dissipative dynamics 
of solitons HHI, in which the "collective coordinate method" is used to transform the 
original Hamiltonian into one of a particle coupled to an infinite set of modes, we show that 
the Bloch wall behaves like a Brownian particle. The advantage of using this method is 
that we keep closer contact with the microscopic details of the system and the mobility is 
naturally calculated as a function of the temperature. The information from the microscopic 
scattering processes between the Bloch wall and the residual modes can be obtained from the 
knowledge of the phase shifts of the associated spectral problem. In the case of reflect ionless 
potentials, as it happens for vanishing anisotropies or external field, the motion of the wall is 



undamped. If this is not the case, the reflection coefficient does not vanish and the mobility 
is finite. 

The outline of this paper is as follows. In Sec. II we present the model. The dynamics 
of its static solution is investigated in Sec. Ill and there we also show how to obtain an 
effective Hamiltonian for the Bloch wall coupled to the residual magnetic excitations. In 
Sec. IV the mobility of the Bloch wall is studied in terms of the scattering phase shifts of 
the second variation problem. The case of large external fields is investigated in Sec. V where 
the phase shifts and the damping constant are explicitly evaluated. Finally, we present our 
conclusions in Sec. VI. 

II. THE MODEL AND ITS STATIC SOLUTIONS 

In this work we consider a one dimensional magnetic system composed by an array of 
spins lying along the ^-direction. Furthermore, let us assume that there is an easy-plane 
anisotropy which tries to keep the spins on the x-y plane and, on top of this, an in-plane 
anisotropy tending to align them along the x-direction. This is a totally anisotropic model 
which is described by a XYZ model of magnetic systems defined by the Hamiltonian 

h = -Y: (J x s\ x) sf ] + JyS^sf + J^sf) -^by: si x) (i) 

(V) h i 

where J x > J y > J z > 0, Sj is the a component (a = x, y, z) of the i th spin of the system, /i 
is the modulus of the magnetic moment on each site and B is the external magnetic field. The 
ferromagnetic XYZ model is actually defined for B = and this is the starting point of our 
analysis. As we can see from (|1|), the ground state of this system is the configuration where 
all the spins are aligned in the ^-direction. However, there is another possible configuration 
which is a local minimum of the energy functional and cannot be obtained from the previous 
uniform configuration by any finite energy operation. 

Let us imagine that we describe our spins classically by vectors 

Sj = S^sin 9i cos ipi, sin 9i sin ifi, cos 9 t ) (2) 



where Qi and ipi are the polar angles of the i th spin. In this representation the above- 
mentioned configuration consists of all Q^s equal to tt/2 and f^s equal to zero or it. 
However, there are other configurations in which Qi = n/2,tfi = if i — ► — oo and fi = tt 
if i — > oo which are approximately (only because Q^s may slightly vary ||) local minima of 
the energy functional of the system. So, Sj winds around the z direction starting at (Q, ip) 
= (tt/2,0) and ending at (Q, <f) = (tt/2,tt). The so-called 7r-Bloch wall § is one example 
of these configurations where ifi varies from to 7r without making a complete turn around 
the z axis. Later on we will see the specific form of this configuration when we consider 
the system in the continuum limit. It will then be shown that Bloch walls are related to 
soliton-like solutions of the non-linear equations which control the spin dynamics in the 
semi-classical approximation. 

If we now turn the external field B on what happens is that the degeneracy between 
if = and if = it is broken. For B > it is clear from (|]) that if = has lower energy 
than if = tt which is now a metastable configuration of the system. In this circumstance 
the system still presents a local minimum of the energy functional. The only difference is 
that whereas Qi is still approximately tt/2, fi starts and ends at zero as -oo < i < oo. The 
27r-Bloch wall is now the configuration where fi winds only once around z. 

In any of the two cases mentioned above, there is no way we could spend a finite energy to 
transform the Bloch wall into the uniform configuration. We would need to turn an infinite 
number of spins over an anisotropy energy barrier. We say that these two configurations are 
topologically distinct. 

Another important point is that ([!]) is translation invariant and this is reflected by the 
translation invariance of the Bloch wall. This means that the region about which the spins 
wind up can be centered anywhere on the z axis. In reality they can even move with constant 
speed along that direction. 

These structures can be obtained by mapping the original Hamiltonian (Jf]) into a 1 + 1 
field theoretical model such as the if 4 , sine-Gordon or any other appropriate model. This 
can be done by simply proceeding a bit further with the semi-classical description for the 



spins we have seen in @ and writing the Hamilton equations of motion for 9i(t) and ifi(t). 
After having done that we take the continuum limit d; t — ► 9(z, t) and <fi(t) — > (f(z,t) and 
write 

e{z,t)n~ + a{z,t), (3) 

where a(z,t) <C 1. Assuming that the variations of if and a from site to site of the spin 
chain are small and linearizing the equations of motion with respect to a one obtains 

if = a2S(J x cos 2 ip + J y sin 2 ip — J z ) (4) 



f)2 

a = a 2 S(J x sin 2 9? + J y cos 2 l f)-^-^ — sin 99 ^- — \- 2S(J X — J y ) cos 99 



(5) 



where a is the lattice spacing. Then, eliminating a from these equations, we get an effective 
equation of motion for (p(z, t) of the form 

1 d 2 ip d 2 ip 



c 2 dt 2 dz 2 



-A\ sin if — A2 sin 2ip (6) 



where 



2a 2 S 2 J x J y ( 1 - ^ ] , (7) 



Jy , 



Ai = -5^T, (8) 



a 2 SJ x h 

and 

^M 1 -!)- (9 » 

Notice that if J x = J y one has A 2 = whereas if -B = it turns out that A\ = 0. So, as 
we can derive the r.h.s. of (H) from a potential energy density U(<f) given by 

U(if) =A 1 (l-cos^) + - i -(l-cos2^) (10) 

we see that Ai controls the potential energy barrier due to the presence of B ^ and A 2 
controls the anisotropy energy barrier. 



The static solutions (dip/dt = 0) of Eq. @ are obtained using that @ 

r (z) d( P' nn 

z-z = / (11) 

are the solitons of the system. In particular, the examples of Bloch walls we gave above are 
the solitons 



ip(z) = 2t&n- 1 [exp^f 2 A 2 (z-z )} if A 1 = (12) 

if (z) = A tan' 1 [exp ^A^z - z )] if A 2 = (13) 

while for the general case of finite anisotropy and magnetic field the solution is the 27r-Bloch 
wall 

i r cosh p , . , , 

<p(z =2tan- x , ; 14 

smh(2/A) 

where we define 



A = I/a/A + 2A 2 (15) 



coshp=^l + ^ (16) 



The soliton (13) can be expressed as a superposition of two twisted 7r-Bloch walls [10 



with arguments \~ l (z — z ) ± p. 

We mention here there is another static solution of (||), the so- called nucleus [|TT| which 
correspond to a superposition of two untwisted 7r-Bloch walls [[TI| . This solution is topolog- 
ically distinct from the previous one. 

III. DYNAMICS OF BLOCH WALLS 

The quantum dynamics of our spin system can be analyzed by studying the quantum 
mechanics of the field theory described by the action 



The next step is to quantize the system described by (|17|). The standard way to carry 
this program forward is to evaluate J7| 

G{t) = tr f Vip exp ~S[<p] (18) 

where the functional integral has the same initial and final configurations and tr means to 
evaluate it over all such configurations. 

As the functional integral in (pl|) is impossible to be evaluated for a potential energy 
density as in ([!(]) we must choose an approximation to do it. Since we are already consid- 
ering large spins (S ^> h/2), and consequently in the semi-classical limit, let us take this 
approximation as the appropriate one for our case. 

The semi-classical limit (h — > 0) turns out to be very easily tractable within the functional 
integral formulation of quantum mechanics J7|. It is simply the stationary phase method 
applied to (0). Moreover, since we are only interested in static solutions, the functional 
derivative of S happens to be the equation of motion (|(J) when dcp/dt = 0. Its solutions can 
be either constant (uniform magnetization) or the solitons (Bloch walls) we mentioned in 
( |T2| - |bl| ). Since we are interested in studying the magnetic system in the presence of walls it 
is obvious that we must pick up one of those localized solutions as the stationary "point" in 
the configuration space and the second functional derivative of ([lOD should be evaluated at 
this configuration. 

When this is done we are left with an eigenvalue problem that reads 

i ah + U " M } ^ z ~ Zo) = ^^ ~ Zq) (19) 

where ip s is denoting the soliton-like solution about which we are expanding ip(z,t). 

Now one can easily show that d(p s /dz is a solution of (|1|) with n n = 0. The existence of 
this mode is related to the translation invariance of the Lagrangian in ( |TTD and this makes 
the functional integral in ( |T8| ) blow up in the semi-classical limit (Gaussian approximation). 

The way out of this problem is the so-called collective coordinate method which was 
developed by field-theorists in the seventies (see |7[] and references therein). It consists of 
keeping the expansion of the field configurations about <f s (z) as 



<f(z, t) =(p s (z- Zq) + Yl Cn^Pn {z - Z ) (20) 

71=1 

but regarding the c- number z as a position operator. Eq. (11) is then substituted in the 
Hamiltonian 

H=J ^ a jJif + u^y +m \, (21) 



where II = -(-£), which can be transformed into [12 1 



H = -^ (P - £ hg mn b + n b m ? + E «Ufr. (22) 

where Q n = ck„. 

In the Hamiltonian fl22|) , P stands for the momentum canonically conjugated to zq, 

M s = —^ / dzU(<p s (z)) (23) 

C J — oo 

is the soliton mass [0 and the coupling constants g mn are given by 



1 

9mi 






dzi} m {z)— . (24) 



2ia 

The operators h + and 6 are respectively creation and annihilation operators for the 
excitations of the magnetic systems (magnons) in the presence of the wall. 

It should also be stressed that Eq. (|22| ) is not an exact result. It is only valid in the limit 
h — ► or, to be more precise, when g 2 h — > where g 2 = -j^ is the coupling constant that 
originally appears in U(g,ip). It must also be emphasized that we have neglected inelastic 
terms such as b + b + or bb because these are only important if the wall moves at high speed 
(v > c) originating Cherenkov-like radiation of the elementary excitations of the medium. 
This approximation also means that the number of excitations in the medium is conserved. 

IV. MOBILITY OF THE BLOCH WALL 

At this point we are ready to start to study properties such as the mobility of the wall 
because we have been able to map that problem into the Hamiltonian fl22|) which on its turn 

8 



has been recently used to study the mobility of polarons, heavy particles and solitons in 
general. We shall not discuss this specific problem in this paper and urge those interested 
in the details of this calculation to follow them in references P,|Il2|-|i"5[ . 



The result that can be obtained reads |II4 



Yl /"OO /"OO 

7(0 = 7777 / / duduj'S(u>, u')(u! — u')[n(u) — n(u')] cos(o; — u')t (25) 

2M Jo Jo 

where 7(t) is the damping function (the inverse of the mobility), 

»(<") = i^n ( 26 ) 

is the Bose function and 

S(u, u)') = J2 \9mn\ 2 5(uj - n n )5(uj' - tt m ) (27) 

ran 

is the so-called scattering function. 

In the long time limit j(t) can, to a good approximation, be written as 

7 (t) = l(T)5(t) (28) 

and j(T) is given by [|I?| 

1 r°° 3E e pE 

W = ^m,1 dEK ^w^w (29) 

where 7Z(E) is the reflection coefficient of the "potential" U"(<p s ) in the Schrodinger-like 
equation ([19|) . Notice that (|29|) is only valid if the states involved in (27) are scattering 



states (see section V.A below for details). One important point that should be emphasized 
here is that there are parameters of the non-linear field equations for which the localized 
solutions render U"(<p s ) a reflectionless potential. These are genuine solitons and for these 
the mobility is infinite. One may realize this is what happens for the Bloch walls ( |I2"D and 
(O). In these cases, the "potential" appearing in flT9|) can be written as 



U"(z) = ?7 2 (1 — 2 sech 2 jjz), (30) 



where f] 2 = A\ for vanishing anisotropy and rj 2 = 2A2 for vanishing external field. The 
spectrum of (p0|) contains a bound state with zero energy 



V>o = J 2 sech(i]z), ^0 = ° 



(31) 



which constitutes the translation mode of the domain wall (Goldstone mode), and a contin- 
uum of quasiparticles modes (magnons) given |16J by 

k n + it] tanh^z) 



^n{x) 



k n +17] 



„ r/Crr >C 



(32) 



where 



2wr S(k n ) 



5(h) = arctan 



2r]k 
k 2 — rj 2 



(33) 



It is known that the reflection coefficient TZ for a general symmetric potential can be 



expressed in terms of the corresponding phase shifts as [J I 



TZ(k) = sin 2 (5 e (k) - 5°(k)) , 



(34) 



where 5 e and 5° are the even and odd scattering phase shifts, respectively. Then, re- 
expressing ( P^ ) in terms of even and odd defined parity states, it is easy to prove that 
this potential belongs to the class of reflect ionless because its phase shifts are given by 



6 e '°(k) = arctan(ry/A;), 



(35) 



that do not distinguish between odd and even parities. 

Nevertheless, when both the anisotropy and external field are finite, the reflection coeffi- 
cient is nonvanishing and consequently the 27r-Bloch wall ([14]) has a finite mobility. In this 
case, the spectral problem ([19]) can be rewritten as 

r d 2 ^ 



dz 2 



+ V(z) >^n(z) = K n ^ n (z), 



(36) 



where the potential V(z) is expressed as 
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V(z) 



sech 2 (p) 

A 1 
tanh 2 (p) 

X 2 



tanh(— + p) tanh(— — p) — sech( — + p) sech( — — p) 
A A A A 

( tanh(— + p) tanh(— — p) — sech( — + p) sech( — — p) 
V A A A A 



y X f Z Z 

sech 2 ( — + p) sech 2 (— — p) I sinh( — + p) sinh(— — p) 
A A V A A 



(37) 



Now, for all finite values of A and p, the translational invariance of the system persists 
and as a consequence the potential (|37D has a zero energy state that is given by 



z z 

ipo oc sech(- + p) + sech(- - p) 
A A 



(38) 



which is nothing but the Goldstone mode of the Bloch walls for finite anisotropy and external 
field. 

In order to obtain an expression for the damping constant (|29|) we need an expression 
for the odd and even phase shifts of (|37D. Unfortunately, their analytical evaluation is very 
complicated for all finite values of A and p, and in what follows we study the situation of 
large fields. 

V. 2tt-BLOCH WALLS FOR LARGE FIELDS 

In this section we evaluate the scattering phase shifts in the situation of large external 
fields and provide an explicit expression for the damping constant. 

A. Scattering phase shifts 



In the case of large fields (p ^ 1) the Shrodinger-like equation (^) can be written as 

r <f 



dz 2 



+ V{z)\i, n {z) = {K 2 n -r 1 2 -P-)i, n {z\ 



t 

A 2 ' 



(39) 



where the potential (|37|) is now reduced to the sum of the reflectionless contribution and a 
perturbation coming from the presence of the large field , explicitly 



V(z) = V (z)+p 2 V 1 (z), 



(40) 
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where 



V (z) = -2sech 2 



(41) 



and 



VUz) =- 1 - 8tanh 2 f ? ] sech 2 ( ? ) . 



A 



A 



(42) 



In order to obtain the even and odd scattering phase shifts A e, ° corresponding to a 
particle in a one dimensional symmetric potential like (f40f ), we will use of a ID version of 
the Fredholm theory |17|], which states that 

,A e '°{E') 



7rA e ' (£)cot(A e '°) 



P / d£- 
h E-E' 



(43) 



where A e, ° are the phase shifts originated by both contributions, the first coming from the 
reflect ionless potential Vb, and the other associated to the high field perturbation V\. On 
the other hand the even and odd spectral functions, A e, °(E), can be calculated from the 
series expansion (see |T7| for details) 



A(E) = -(E\V(z)\E)+V f 

Jo 



°° dE 1 



E — E\ 



(E\V(z)\E) {E\V(z)\Eh) 
(E^V^IE) {E x \V{z)\E x ) 



+ 



(44) 



where V stands for the Cauchy principal value. Clearly the expression (|44f) cannot be 
analytically evaluated to all orders. On the other hand, making use of (|35|) and considering 
that p is small enough, the expression for the the phase shifts (|43"D can be written up to first 
order in p 2 as 



V , 2 ^ A i 



^'■° = i + r\ + 2Br 



(45) 



where 



A 1 = -(E\V 1 (z)\E), 



(46) 



and 
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5n 



V 



- A*'°{k')k'dk' 

-co k 2 - k' 2 ' 



A) = -(s|y (^)|s). 



Using a convenient basis set, the three expressions given by 
lytically evaluated (see the apendix) yielding 



and 



A e '°(k) 



2M 



± 



nkX sinh(7rA;A) 



AT{k) 



8p 2 M 



' T(2k 2 -\- 2 ) ' J 



nk\ 



sinh(7rA;A) 



B e,o 



AM °° 

— E(-i) n+1 



n 



h 



71=1 



((k/X) 2 + n 2 



(47) 
can be ana- 

(48) 
(49) 
(50) 



Now we can finally write down an expression for the phase shifts A e, ° by substituting p9] 
and (|5T)D in (^5|). In so doing one gets 



tanA e '°(A;) = ? + 



-2 



87rh~ z Mp 



i i 8MV/ T>n+1 w 



T 



2A; 2 A 



2\2 



A;7rA 3sinh(A;7rA) 



(51) 



(k/X) 2 + n 2 

Looking at (51) we realize that, whereas A e (k) remains almost unchanged as a function 
of k for p / (see Fig.l), A°(k) presents a completely different structure. As it is shown 
in Fig.2, the odd phase shift for p ^ display a rr amplitude discontinuity characteristic 
of a resonance. This discontinuity becomes sharper as p —>■ 0, giving no contribution to 
the reflection coefficient, because as expected from (|55]), the even and the odd phase shifts 
approach each other for small p (see Fig. 3). In this situation the ir amplitude discontinuity 
is still present, but as an isolated one, does not contribute to the reflection coefficient. We 
can now turn our atention to the behaviour of the phase shifts near k = 0. In order to do 



this we will use the ID version of the Levinson's theorem |[J| for one-dimensional symmetric 
potentials. Levinson's theorem establishes that: 

1 



A e (k = 0) = n(n e - 
A°(k = 0) = im°, 



(52) 
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where n e and n° are the number of even and odd parity bound states. As it can be seen in 
Fig. 1, the even parity state phase shift is always 7r/2 for k = for any finite value of p, 
then from ([52]) it is clear that our system has an even parity bound state for any value of 
p. The existence of this state is in complete agreement with the translational invariance of 
the system and corresponds to the Goldstone mode. The analysis of the existence of an odd 
parity bound state is slightly different. As can be seen for (^) the correct definition of 5° 
according to (|||) is 



6°= I 



(53) 



0, k = 

arctan(?7//c), k ^ 

because there is no odd parity bound state for the p = case. Following the same idea 
and realizing that ([51]) behaves in the same way as (|35|) for small enough k, we conclude 
that A°(0) = and therefore there is no odd bound state in the p ^ case. Actually, 
A e (0) ~ A°(0) for large values of k whereas they start to deviate from one another as k 
decreases. Equation (^3|) represents an extreme situation when 8°{k) discontinously jumps 
from 5 e (k) to 5°(k) = when k —» 0. This is the only way we can reconcile the absence of 
an odd bound state and the reflectiolessness of the potential. Therefore, the spectrum of 
(pfl]|) is composed by: i) the ip solution (^) corresponding to the translation mode of the 
wall (Goldstone mode) and ii) the ipk solutions which constitute the continuum modes and 
correspond to magnons. 

B. The damping coefficient 

In order to find the damping coefficient we must compute lZ(k). This can be done by 
inserting ([51]) into the general expression 

K{k) = sm 2 (A e (k) - A°(k)). (54) 

Although we have its analytical form in the present approximation, we had better plot it for 
the whole range of the momentum k for various ratios of A 2 /Ai (anisotropy /external field) 

14 



(see Fig. 4). Having done that, one can immediately integrate this function in expression 
(p9|) which finally allows us to describe the damping as a function of the temperature as 
shown in Fig. 5. As it can be seen, the damping constant is linear for the whole range 
of temperatures. This result can be obtained directly from fl29|). For T high enough the 
damping constant can be approximated by 

which is linear on T, independently of the explicit form of TZ(E). In the low temperature 
regime we can write 

1 POO 

^2^/„ USnWB*-*' (56) 

where E always presents a gap given by Eg = [(t]X) 2 + p 2 ]h c 2 /A 2 and therefore and expo- 
nentially small damping. However, this expression must be studied in the limit in which 
we are interested; namely, p — > 0. Unlike 7(T) given by ([55]) the expression (|56| ) cannot be 
analitically estimated in a trivial way because TZ{E) has only been computed numerically. 
Here it must be stressed that this result is only reliable for not too low temperatures because 
we have employed the odd phase shift A°(k) from Fig. 2 and its computation clearly do not 
account for its correct values for k — > 0. This behaviour is due to the fact that for very low 
energies one can not approximate A(E) in (|44]) by its first term although this approximation 
works well for the even phase shift A e (k). This turns out to be a good description of the 
mobility of the wall for extremely high fields if we keep the above explanation in mind. 

VI. CONCLUSIONS 

In this paper we have analyzed the dissipative dynamics of Bloch walls in a one dimen- 
sional anisotropic ferromagnet in the presence of an external magnetic field. In particular, 
we have considered the limit of high magnetic fields although there is no reason why one 
should not apply the same methods to the low field case. The only difference is that the 



15 



scattering problem with which one has to deal is more straightforward in the high field 
case. Our predictions are that the damping coefficient j(T) presents a linear behaviour as 
is plotted in Fig. 5, for a vast range of temperatures. 
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VIII. APPENDIX 

In this appendix we show how to obtain the expressions of the even and odd spectral 
functions ( f47|) and (f46|) . Suppose we have a particle in a one dimensional symmetric potential 
of the form V = Vq + gV\ confined to a region (— L, +L) with L much larger than the range 
of the potential V. The asymptotic form of the wave functions for V 7^ are given by 

\z) e = JjCos(k\z\+A e (k)), 

\z)° = J \ sgn(z) sin(*|z| + A°(*)), (57) 

for \z\ — > oo. If V = the wave functions have the same structure as in fl57|) with A e, ° = 0. 
Because the wave functions must vanish at z = ±L one realizes that 

SE n _l_ Ae , 



AE n n 

where AE n = E® +1 — E® and 5E n = E n — E®. Following closely the prescriptions given in 



17| for the 1-D case, the spectral functions (j47|) and (|4q ) are given by 

16 



/+oo 
dzV t (z)(z\E) e , (E\z) e , , 
-oo 

where the states (E\z) £)0 can be obtained by taking 



i = 0,l 



(z\E) et0 = lim 



Explicitly, 



(z\E) 



L^OC y/AE n 



cos(kx), for even parity 



y Ltxk sin(foc), for odd parity 
Inserting (|60|) and ([42]) in (|58"D we have 

A? + A? = ^^ / + °° sech 2 (|) tanh 2 (|)dz, 

which can be easily evaluated with the substitution y = tanhz/X, yielding 

16p 2 M 



A\ + A\ 

On the other hand we have 

8Mp 2 r+oc 



nh 2 \k 



A e — A° 

7m K J-oo 



sech 2 (— ) tanh 2 (— ) cos(2k\)dz, 
-oo A A 



that can be writen as 



A\- A\ 



o 8Mp 



2 ,+ 



Trr/c 



/+oo f z Z \ 

sech 2 (-) - sech 4 (-) cos{2kX)dz 
-oo V A A / 



which can be analytically evaluated [[H| yielding 

16p 2 M 



A^-A° 



h sinh(7r/cA) 



1 2fc 2 A 2 
3 3 



(58) 



(59) 



(60) 



(61) 



(62) 



(63) 



(64) 



(65) 



Therefore combining (|5"2]) and (BH) we have PO). In the same fashion it can be shown 



that 



A° + A° 



AM 



and A e - A° 



AM 



(66) 



which immediately gives ([|8 ) . Now we can evaluate the Cauchy principal value in Q4"7| ) which 



reads 



17 



r,eo 2Af r+^ 



1 1 

± 



irk'X sinh(7r/c / A) 



k'dk' 
k 2 - k' 2 ' 



(67) 



The first term on the right hand side of fl6*7|) is clearly zero. Therefore, using the product 



expansion of sinh(7rz) fuction [19] its second term becomes 

dq 



R e,o 



2M '■+--■ 



n 2 



n 



(68) 



nh 2 ' J-oo k 2 — q 2 ^ n 2 + q 2 

where q = k/\. Going to the complex plane, the previous expression can be analitically 
evaluated as 



R e,o 



4M ™ 



-l) n n 



V ^i (k/X) 2 + n 2 



(69) 
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FIGURES 



FIG. 1. The even phase shift as a function of the momentum for three different situations. 
The continuous line corresponds to p = 0.14, the dotted line to p = 0.31 and the dashed line to 
p = 0.60. 



FIG. 2. The odd phase shift as a function of the momentum. The triangles correspond to 
p = 0.14 and the circles to p = 0.31. 
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FIG. 3. The even and odd phase shifts when p = 0.14. As can be seen they approch each other 
as the ratio A\j A<i increases, thats means p — ► 0, and the only difference conies from the singular 
point in the odd phase shift contribution. 



FIG. 4. The reflection coefficient as a function of the momentum. The continuous line for 
p = 0.14, the dotted line for p = 0.31 and the dashed line for p = 0.60. 



22 



FIG. 5. The damping coefficient as a function of the temperature for different values of p. The 
continuous line for p = 0.14, the dashed line for p = 0.31 and the dotted line for p = 0.60. 
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